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Abstract 


The present work follows a recent survey of airframe noise prediction methodologies. In that survey, Lighthill’s acoustic 
analogy was identified as the most prominent analytical basis for current approaches to airframe noise research. Within this 
approach, a problem is typically modeled with the Ffowcs Williams and Hawkings (FW-H) equation, for which a geometry- 
independent solution is obtained by means of the use of the free-space Green function (FSGF). Nonetheless, the aeroacoustic 
literature would suggest some interest in the use of “tailored” or exact Green’s function (EGF) for aerodynamic noise prob- 
lems involving solid boundaries, in particular, for trailing edge (TE) noise. A study of possible applications of EGF for 
prediction of broadband noise from turbulent flow over an airfoil surface and the TE is, therefore, the primary topic of the 
present work. Typically, the applications of EGF in the literature have been limited to TE noise prediction at low Mach 
numbers assuming that the normal derivative of the pressure vanishes on the airfoil surface. To extend the application of EGF 
to higher Mach numbers, the uniqueness of the solution of the wave equation when either the Dirichlet or the Neumann 
boundary condition (BC) is specified on a deformable surface in motion. The solution of Lighthill’s equation with either the 
Dirichlet or the Neumann BC is given for such a surface using EGFs. These solutions involve both surface and volume 
integrals just like the solution of FW-H equation using FSGF. Insight drawn from this analysis is evoked to discuss the 
potential application of EGF to broadband noise prediction. It appears that the use of a EGF offers distinct advantages for 
predicting TE noise of an airfoil when the normal pressure gradient vanishes on the airfoil surface. It is argued that such an 
approach may also apply to an airfoil in motion. However, for the prediction of broadband noise not directly associated with 
a trailing edge, the use of EGF does not appear to offer any advantages over the use of FSGF at the present stage of develop- 
ment. It is suggested here that the applications of EGF for airframe noise analysis be continued. As an example pertinent to 
airframe noise prediction, the Fast Scattering Code of NASA Langley is utilized to obtain the EGF numerically on the 
surface of a three dimensional wing with a flap and leading edge slat in uniform rectilinear motion. The interpretation and 
use of these numerical Green functions are then discussed. 
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1. Introduction 


In the present work we continue the study of airframe noise (AFN) prediction methodologies ^ . This investigation is moti- 
vated by the need to develop physics-based AFN prediction tools. As future airport-noise abatement goals become more 
stringent, it will eventually become imperative to predict overall noise levels during an aircraft’s design stage, including its 
airframe design. The success of such a “quiet-design” approach will depend on a physics-based prediction methodology, i.e., 
a methodology in which the noise prediction process incorporates modeling that is derived from first principles. Noise 
prediction tools that include as much physics-based modeling as possible are highly tractable and robust, and thereby well 
suited for future design guidance. 

The current AFN methodology study began with a survey of the prediction methods that are currently in use or in various 
stages of development 1 . In that work, present approaches to AFN prediction are categorized into four groups, three of which 
are referred to as “semi-empirical,” “fully analytic,” and “fully numerical.” In the fourth approach, a computational fluid 
dynamics (CFD) simulation provides the source description, and the radiated sound is predicted by an appropriate acoustic 
formulation based on Lighthill’s acoustic analogy . At present, only the CFD/acoustic- analogy (CFD/AA) approach pro- 
vides the potential to develop the necessary physics-based prediction tools within the required time frame ^ . As this conclu- 
sion has remained relevant for the past decade, NASA’s development of a physics-based AFN methodology has been largely 
concentrated within the CFD/AA approach. 

The advantages of a CFD/AA methodology result from the very notion of “acoustic analogy,” itself. Morris and Farassat 
describe an acoustic analogy as “any aerodynamic noise theory in which the equations of motion for a compressible fluid are 
rearranged in a way that separates linear acoustic propagation effects.” ^ Furthermore, “this rearrangement results in a set of 
equivalent sources that are assumed to be non-negligible in a limited region of space.” Therefore, within the CFD/AA 
approach, the CFD solution is regarded as input data to a stand-alone acoustic propagation algorithm. As such, the propa- 
gated sound can be predicted to observers in arbitrary locations, e.g., at distances many times larger than the spatial domain 
of the input CFD solution. Moreover, such an acoustic-analogy based algorithm is a direct application of an explicit, analyti- 
cal formulation for radiated sound. Thus, there is no additional error incurred during the acoustic calculation; i.e., the far- 
field noise prediction is as accurate as the computed source input. 

The most commonly used governing equation for AFN prediction is the Ffowcs Williams-Hawkings (FW-H) equation^. The 
solution of this equation is obtained by using the free-space Green’s function (FSGF) which is independent of the aircraft 
geometry. The usefulness of this approach has been clearly established by the prevalence of noise prediction tools based on 
solution of FW-H equation. Nonetheless, in recent years a number of technical papers and books have appeared that promote 

C_ 1 o 

the use of the so-called tailored or exact Green’s function as an analytical tool for AFN problems . Although Howe has 
been the leading advocate of the use of EGF, the present authors believe that it was Goldstein who first suggested the use of 
EGF in an important paper in 1974^ and in his book on aeroacoustics ^ . In Goldstein’s works, the EGF was discussed in 
connection with the possible choices of the Green’s function in acoustic problems that are solved in the presence of solid 
bodies. 

We feel that it is time that the use of EGF approach be reviewed carefully and compared to the use of the FSGF in AFN 
prediction problems. Such a review and comparison are presented in this paper. Our approach is as follows. In Section 2, we 
briefly discuss the history of broadband noise prediction problem, what is known from previous research, and why the use of 
EGF has attracted so much attention from aeroacousticians. Then we explicitly state the mathematical problem that we 
address here which is whether EGF approach can be used for a surface in motion. We accomplish this by proving some 
uniqueness theorems for the wave equation. To this end, in Section 3 we prove rigorously that the solution of the initial- 
boundary value problem of the wave equation is unique with either the Dirichlet or the Neumann boundary condition speci- 
fied on a deformable body in motion. This means that the EGF may be found for a moving surface and used for AFN predic- 
tion. In Section 4, we discuss the implications of this result when used for broadband noise prediction. Our main conclusion 
is that there seems to be clear advantages of using EGF over the FSGF for TE noise prediction at low Mach numbers and 
perhaps even at high Mach numbers. It is recommended that the use of EGF be explored further in AFN prediction. Finally, 
in Section 5 we give examples of EGF for the nontrivial case of two three dimensional wings one of which has a trailing edge 
flap and a leading edge slat. These EGF are constructed numerically for a fixed source frequency on the wing surface for 
several observer positions . Concluding remarks follow in Section 6. 
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2. The Airframe Broadband Noise Problem 


2.1. Some History 

Broadband noise radiation from airfoil surfaces has been studied in connection with aircraft and underwater acoustics by 
Amiet and coworkers , Howe ’ , and Brooks and coworkers since the early 1970s. Until recent years, 

almost all of noise prediction methodologies depended on the use of Ffowcs Williams and Hawkings (FW-H) equation ^ 
whose solution uses the FSGF of linear wave equation. The use of FSGF allows various closed form solutions of FW-H 
equation because of extreme simplicity of the analytic form of FSGF and its easily visualizable geometry, e.g., the collapsing 
sphere visualization of the characteristic cone of the wave equation in four dimensions. Two fundamental papers by Ffowcs 
Williams and Hall^ and Crighton and Lippington^ have contributed greatly in understanding the physics of TE noise 
generation by quadrupoles and vortices passing in the vicinity of the edge. The results of these authors have been verified by 
many others . 

In recent years researchers within NASA, the universities, aircraft and engine industry in the U.S., as well as similar organiza- 
tions in Europe and Asia have conducted intensive research on AFN. Many technical papers have been published. There has 
been a sudden surge of interest in the use of EGF as evidenced by the number of authors who are using EGF in their publica- 
tions. The chief advocate of the use of EGF in recent years is Howe^’^ . But Goldstein had already suggested the use of EGF 
in a paper in 1974^ and in his well-known book aeroacoustics ^ which is based on an earlier NASA publication with the 
same name^ . It is also true that the authors of the two fundamental papers on TE noise^’^ used the EGF for the two 
dimensional problems that they considered. We will discuss below some of the advantages of the use of EGF in AFN 
prediction. 

An open and important question that the researchers who propose the use of EGF have not addressed involves the uniqueness 
of the solution of the wave equation with the Dirichlet or Neumann boundary condition (BC) on a moving surface. Appar- 
ently, the prevailing view is that because of the slow speed of aircraft during take-off and landing (about Mach 0.2), the 
uniqueness theorem for a stationary surface will apply in AFN predictions. Although, intuitively, the solution’ s uniqueness 
appears obvious, it is not enough to rely on intuition alone. If the uniqueness theorem breaks down for a moving surface, 
much of the recent results based on EGF will be in doubt. We, therefore, propose to prove here two uniqueness theorems for 
the wave equation with the Dirichlet or the Neumann BC for a moving surface. These two theorems constitute the main 
contributions of the present paper. We will then discuss the implications of these results as used in broadband noise predic- 
tion. 


2.2. What We Know Now 


We give a summary of what is known in predicting broadband noise caused by turbulent flow over an airfoil and the TE 
noise. Let us start with the acoustic analogy using the FW-H equation whose solution requires FSGF. The FW-H equation as 
originally proposed by Ffowcs William and Hawkings is^ : 


8 


d 


□ P ' = ~x~ [Po 6(f)] - — — [p rii 6(f ) ] + 


dt 


dx; 


d 2 

dxj dx. 


-[ H(f)Ty] 


( 1 ) 


Here, the moving surface is described by f(x , t) = 0 such that vf = n , n is the unit outward normal. This assumption implies 
that / > 0 outside the moving surface (see fig. 1). Also p' = c 2 p' = c 2 (p - p 0 ) , c and p 0 are the speed of sound and density 
in the undisturbed medium, respectively, v n , p and Ty are the local normal velocity of the surface, the local gage pressure of 
the surface, and the Lighthill stress tensor, respectively. The Heaviside and the Dirac delta functions are denoted H(f) and 
6(f) , respectively. 

Ffowcs Williams later proposed to use a penetrable (porous or permeable) data surface to account for nonlinearities in the 

7Q 

vicinity of a moving surface^ . We again assume that the penetrable surface defined by / (x, t) = 0 and the fluid velocity is 
denoted by u. The FW-H equation for penetrable data surface, FW - H v ds , is^ y ^ : 
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□ 2 cl p' = n 2 p' = A[ Po u n 6(f)] - ~[U 6(f)] + -AF- [Tij H(f)] 
We have used the following notations in the above equation: 



( 2 ) 

( 3 ) 


U = p Sij rij + p Ui (u n - v n ) (4) 

where S t j is the Kronecker delta. The philosophy behind using FW-Hp^ is to locate the data surface / = 0 to enclose a 
moving noise surface, in such a way that all quadrupoles producing non-negligible noise are included within this surface, so 
that no volume integration of the quadrupoles outside the data surface is necessary. This data surface must be located in the 
region of high grid density in order to adequately account for the upper spectrum of a broadband AFN calculation. For this 
reason, one would like this surface to be as small as possible, because of the computer intensive nature of turbulence 
simulation^ l - 36 

The formal solution of eqs.(l) and (2) has two surface integrals involving local normal velocity and surface pressure and a 
volume integral of the Lighthill quadrupole^’30 . This solution form is a consequence of the fact that we have used FSGF in 
the derivation. For studying broadband noise some sort of statistical analysis of the solution must be performed as, for 
example, is done by Morris and Farassat for jet noise analysis ^ to relate the statistics of the noise to the statistics of turbulent 
flow. When both surface and volume integrals appear in the solution, the algebra for statistical analysis becomes very 
complicated. 

Concerning the surface data needed for noise calculations, it is clear that, because of the appearance of the Dirac delta 
functions in eqs. (1) and (2), no normal derivative of surface pressure should appear in the solution. This is expected because 
of the following result from generalized function theory-^ ' ’■ : 

p (x, t) 6(f) = [p (x, t)] f=0 6(f) => v[p (x, t)] f=0 = v r p (x, t) (5) 

where V T p (x, t) is the gradient of p (x, t) tangent to the surface. Therefore, the surface input data used for source descrip- 
tion for FW-H equation is of Dirichlet type. It is important to remember this when we discuss the uniqueness theorems for 
the solution of wave equation with the Dirichlet or the Neumann BC. 

Some important fact to know about the use of FW-H and equations for broadband noise prediction are: 

1- These equations are exact and given the exact input data that accurately models the source, one should get the correct 
radiated noise , 

2- The surface pressure p may include an acoustic component which should be included in predicting noise off the surfa- 
4 97 

ce^ 4 * * 7 (see next subsection), 

3- There is the possibility of noise cancellation from surface and volume sources in eqs.(l) . This phenomenon is direction 
and frequency dependent. 

4- Doing statistical analysis on solutions based on FW-H and FW - H p ^ s is difficult if all the integral terms in the solution 
are retained. Relating the statistics of the broadband noise to the statistics of the volume and surface sources is essential in 

understanding and reduction of AFN. 

It is important to recognize that, in general, the volume sources involved in the generation of broadband noise occupy small 
region of space in the vicinity of the surface but can have very high velocity gradients. Conceptually, a compressible turbu- 
lence simulation would model the acoustic waves propagating over the surface. The problem is that compressible turbulence 
simulation is computer intensive and thus expensive and we prefer incompressible turbulence simulation. However, at 

present, compressible turbulence simulations at meaningful scales are prohibitively expensive, and therefore research inter- 
ests are pursued in the area of incompressible simulations. The question is whether surface integration of the acoustic pres- 
sure can be avoided if at all possible. The affirmative answer for some noise prediction problems is given in the next 
Subsection. 
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The possible cancellation between volume and surface integrals is another matter. In general, this problem should not exist if 
the penetrable data surface for eq. (2) encloses the entire region of non-negligible turbulence intensity and vortical flow. This 
is not always possible because of the high cost of turbulent flow simulation. Although a limited forms of solutions of FW-H 
equation employing the FSGF have been tested so far, the possibility of adding some of the surface terms hidden in the 
volume source term of eqs.(l) and (2) to the loading term as discussed by Farassat, Myers and Brentner^’^l to avoid this 
problem has not been explored. Also, improvements in flow simulation to extend the region of computation for use in a 
solution based on FSGF in FW - H pds will alleviate the cancellation problem. However, the issue of specifying statistics on 
the data surface, based on the flow simulation, remains. The use of EGF can avoid this difficulty in statistical formulation in 
certain situations, e.g., TE noise prediction, but not the need for an accurate turbulence simulation. 


2.3. Why Use EGF? 

Howe^ discusses the Curie problem^ and the solution that Powell gives using the image of quadrupoles reflected in the 
plane boundary^. Howe points out that to predict the noise using the Curie formula, one has to integrate the acoustic 
pressure over a large portion of the plane boundary. In the Powell formulation, this surface integral is simply equal to the 
contribution of the image quadrupoles. Howe considers Powell’s solution which is based on the EGF more appropriate. We 
will give another example which is the piston in the wall problem. 

Figure 2 shows a piston in an infinite baffle with normal velocity distribution v n (x 1? jc 2 , t) on the piston surface fl 5 . The 
solutions of this problem by using FSGF in FW-H equation, and by Rayleigh using two different EGF based on the Dirichlet 
and Neumann BCs are^ : 


f Po v n r 

4tt/7'(x, t)= - dS + 

Jn , r Jr 


[Het C ° S0 


dS + 


/ 

Jr 2 


[p] ret cos 6 


dS (FW - H eq.) 


2 7T p 


2 71 p 


'(*,*) = f 

Jr 

'(*,*)= f 


[Het C 0 S C„. f Het C0S « 


dS + 


2 cr 


[pov„] 


f 


2 


dS- 


i 


[dp/dn\ 


dS (Dirichlet BC) 


dS (Neumann BC) 


( 6 ) 

(7) 

( 8 ) 


It is clear that the solutions given by eqs.(6) and (7) which require the knowledge of the surface pressure over the entire two 
dimensional space, have basically two shortcomings as compared to that given by eq.(8). First more input data are needed to 
get the acoustic pressure and, second, the prediction is more costly. Note that the simplest solution, eq.(8), uses the Neumann 
BC. 

Let us now discuss why Howe recommends the use of the EGF with the Neumann BC dp/ dn-0 for the TE noise problem 
when the airfoil moves at low Mach number ^ . Let the typical wavelength of the sound associated with the vortical flow over 
the TE be much greater than the chord. Then one can use incompressible turbulence simulation near the TE and use the 
volume quadrupoles as the source of sound which will be scattered by the TE. The use of the Neumann BC eliminates the 
need for the knowledge of the pressure on the surface. This type of BC means that dp / dn is specified on the airfoil surface 
which for a stationary surface is dp / dn- 0. The condition imposed on the EGF in this case is that dG/ dn - O.The only 
sources remaining are the Lighthill quadrupoles which are volume sources. An additional advantage here is that by having 
only one integral to evaluate the radiated noise, statistical analysis of the acoustic pressure based on the statistics of the 
turbulent flow is considerably simplified. Even for the case of high frequency TE noise where the typical wavelength of 
sound is much smaller than the chord, or when the airfoil is in motion, the assumption of dp /dn = 0 appears to be accept- 
able on physical grounds. The turbulent flow simulation near the TE may have to be compressible, but limited to the TE 
region, and the noncompactness of the volume sources may have to be taken into account for noise calculations. 

To study the use of EGF for surfaces in motion, we will next define the problem that we will address in the next section. 
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2.4. Problem Definition 

The broadband airframe noise problem that we will address here is due to the turbulent flow in the boundary layer (BL) and 
the passage of the vortical motion in the BL over the TE of a wing and its control surfaces. As is well-known, these are not 
the only sources of broadband noise of airframe ^ . The turbulent pressure fluctuations on the wing surface and the scattering 
of the energy in the vortical flow over the TE generate broadband noise that we are concerned with in this paper. Theoreti- 
cally, the quadrupoles in the BL and the wake of the wing and control surfaces as well as in the flow over the side edges of 
flaps and slats also generate broadband noise. We are interested in proposing a suitable formulations to predict this broad- 
band noise assuming that high-resolution turbulent simulation data are available to us. In particular, the use of an EGF 
formulation for this problem will be studied. 


3. Some Uniqueness Theorems for the Initial-Boundary Value Problem of 
the Wave Equation 


The problem that we are considering here is described as follows. Let an arbitrary deformable surface moving in the 
unbounded three dimensional space 1R^ be given implicitly by /(x, t) - 0 such that V / = w, where n is the local outward 
pointing unit normal to the surface (see fig. 3). The symbol Fl t stands for the region exterior to the surface / = 0, and dkl t 
refers to the boundary of Fl t described by the equation f - 0 . The fact that n points into the exterior of / = 0 implies that 
/ > 0 in . We want to find out whether the following problem has a unique solution: 

□ 2 0 = 0 (x, t), (x, t) e£l t x [0, T ] 

IC : 0 (x, 0) = g(x), d t 0 (x, t) = h(x), xe£l t 

<90 

BC : either 0 (x, t) - k (x, t ) or — — = k (x, t ) is specified for (x, t)e dkl t x [0, T] 

on 


In this section, we are taking c - 1 in the wave operator, i.e., U 2 = <9 2 -V 2 . 

As is well-known from the theory of partial differential equations, the above BCs are called the Dirichlet and Neumann 
boundary conditions, respectively. We require that the IC and BC be consistent, e.g., for Dirichlet BC, one should have 
k (x, 0) = g (x) on dkl t . 

To prove uniqueness, we must show that if (f) l (x, t) and 0 2 (x, t) are two solutions of eq.(9), then their difference 
0 — 0 1 -0 2 which satisfies 

□ 2 0 = 0, (x, t)ekl t x [0, T] 

IC : 0 (x, 0) = 0 ,<9,0 (x, 0) = 0, xekl t ^ 

<90 

BC : either 0 (x, t) = 0 or — — = 0 for (x, t)e dkl t x [0, T ] 


has the solution 0 (x, t) = 0 for (x, [0, T ] . 

The general method of proof is based on an energy relation. This is described in Subsection 3.2 below. 


3.1. The Classical Uniqueness Theorem 

The classical uniqueness theorem states that for a stationary surface /(x) = 0, the solution of eq. (1) is unique^ . We will not 
give the proof here because it is a special case of our more general result below. 


3.2. Extension of the Uniqueness Theorem to a Moving Deformable Surface 

We want to show that the solution of problem described by eq.(10) is 0 (x, t) = 0 for (x, t)ekl t x [0, T ] . We start with the 
following easily established identity: 
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<!>, n 2 <P = Yt{\ $) “ • V(f> = A (I + J_ | pj _ v . ( 0 t V 0) 

We now integrate this result for (x, t)e£l t x [0, T] and use the divergence theorem. This gives 

ff §^ff a 

^0 UL ^0 


dcj) 

dn 


dSdt 


( 11 ) 


( 12 ) 


where £ = y 0^ + j- | V0 | z which can be interpreted as an energy quantity. We next obtain the following result to use on the 
left side of the above equation: 


d 

dt 


f E(x, t) dx-^~ f H(f) E(x , t)dx- f 
Jf L dt J R 3 


dE(x, t ) 

H(f) -^T- 1 dx 
dt 


-L 


v„ 6(f) Edx 


(13) 


Here H(f) and 6(f ) are the Heaviside function and the Dirac delta functions, respectively. We have also used the relation 
df / dn = — v„ where v„ is the local normal velocity on the surface / = 0. We next use the following relation proved in 
references ’ in the above equation: 


f v n 8(f)Edx= f ' 
Jr 3 Jdn t 


v„ EdS 


(14) 


Equation (13), with some rearrangement, will then yield: 


X. 


H(f) 


dE(x, t ) 
dt 


dx - 


r 8E dr r 

I — — dx- — I Edx+ ' 
Jn t dt dt J nf J d£lt 


v„ EdS 


(15) 


This is, of course, a transport theorem. We next use this result in eq.(12) and keep the surface integrals on the right of equal- 
ity sign to get: 




dSdt 


(16) 


We next integrate the time integral on the left side. This gives 

r^T 


f Edx- f Edx = T f b,!^-v n A 
dn t (T) dn t (0) do dsn A ° n ) 


\dSdt 


(17) 


Here f \(T), and H r (0) are the regions exterior to the surfaces /(x, T) - 0 and /(x, 0) = 0, respectively. Now we utilize the 
IC of eq.(10) in the second integral on the left side of eq.(17) to get 


/ 


Edx - 0 


(18) 


Equation. (17), therefore, becomes 

f Edx= f f U t ^--y n E\ 

dn t (T) Jo Jdn t V dn ) 


dSdt The Energy Identity 


(19) 


This is the generalization of the energy identity used for obtaining the classical uniqueness theorem^ . We now prove some 
very useful uniqueness theorems based on this identity for moving and deformable surfaces for our applications. 

3.2.1. The Uniqueness Theorem for an Infinitely Thin Deformable Surface in Motion 

For such a body v„ = 0 so that eq.(19) becomes identical to what that is obtained for stationary surfaces^ : 


f Edx= f f 

Jn t (T) Jo Jon t dn 


dSdt 


( 20 ) 
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Now using either the Dirichlet or the Neumann BC of eq.(10), the right side of eq.(20) vanishes and we obtain: 


/ 

Jn t (D 


Edx - 0 


( 21 ) 


Since E is a positive quantity, this relation implies that E - 0 . From this we conclude that 0 = constant = 0(x, 0) = 0 . 


3.2.2. The Uniqueness Theorem for an Arbitrary Deformable Surface in Motion With Dirichlet Boundary 
Condition 


This theorem is considerably more difficult to prove than for an infinitely thin moving body. We will again utilize eq.(19) but 
this time some concepts from differential geometry must be used to prove uniqueness of the solution. We mention here that 
by using eq.(19), it is a trivial matter to show that we have uniqueness if consistent Dirichlet and Neumann BCs are specified 
on a moving surface. In fact, in this case we have a closed form solution which is the Kir chhoff formula for moving surfaces 
first derived by Morgans^ i n 1930 and verified much later using generalized function theory by Farassat and Myers^7>48 
We know that in the classical Kirchhoff formula, the imposition of the Dirichlet and Neumann BCs is an overspecification of 
data (from classical uniqueness theorem given above). We also feel intuitively, based on the fact that there is total destructive 
interference of acoustic signal inside the moving surface, that giving the Dirichlet and Neumann BCs on a moving surface is 
also an overspecification. Now that many researchers are advocating the use of EGF for TE noise prediction, there is a need 
for uniqueness theorems of wave equation for either Dirichlet, Neumann on a moving and deformable surface. We now 
proceed with the proof of the uniqueness theorem when only the Dirichlet BC is specified on the moving surface. 

First we find the expression for the integrand of the surface integral on the right of eq.(19) in terms of the quantities known 
on the surface. We then proceed to show that the integral on the right of this equation vanishes in the vicinity of any time 
t = to . This process then can be repeated for all t e [0, T] . This leads to the conclusion that integrand of the integral on the left 
of eq.(19) is zero. The conclusion of 0 = 0 follows by the same reasoning as used in the case of infinitely thin surface. 

Note that eq.(19) was derived in a frame fixed to the undisturbed medium. This means that in this equation <fi t = <90( x, t)/dt. 
This is not the same time derivative of 0 on the moving surface. Referring to fig. 4, we define a Gaussian coordinate system 
(w 1 , u 2 ) on the surface / = 0. Next define the coordinate u 3 locally normal to this surface with positive direction along the 
direction of n with u 3 as the distance from the surface (time t is frozen). Since n = V/, this means simply that we have 
defined u 3 = /,i.e., dn = df = du 3 . We will use the notation u = (u l , u 2 , u 3 ). We can then show that, if we define 
0 (w, t) = 0(x (w, t), t ) 


f t ( 11 , t) = 


<90 (w, t) 
dt 


<90(x(w , t), t) 
dt 


= 0* (x, t) + v-V0(x, t) — (p t (x, t) + v • V 0 (w, t) 


( 22 ) 


where v = dx(u J , u 2 , 0, t)/dt is the local surface velocity which is unambiguously defined. Note that the Dirichlet BC can 
now be written simply as 0 {u 1 , u 2 , 0, t) = 0. This means that V 0 {u 1 , u 2 , 0, t ) = f n (u 1 , u 2 , 0, t)n = f n (x, t ) n where 
0 n = <90 / du 3 . Equation (22) and the BC then give 

ft (X, t ) |y- 0 = ~V n f n (u 1 , U 2 , 0, t) = ~V n f n (X, t) | /=0 (23) 

Using all these results in the integrand of eq.(19), we get 

{$, - v„ #) = — v„ (3 + V 2 n )f (u 1 , u 2 , 0, t) (24) 


The most important thing to recognize at this point is that the expression on the right of eq.(24) is defined in the vicinity of 
the surface / = 0 at any time as a function of coordinates (w, t ) and not just (w 7 , u 2 , 0, t). We use this fact in eq.(25) below. 
Now we will take two steps that are crucial in the proof: 

1- At the time to e [0, T], set up the coordinate system u on the surface /(x, to) = 0 as described above with (u 1 , u 2 ) the 
Gaussian coordinates on this surface and u 3 coordinate normal to the surface, u 3 - 0 on the surface and positive in the 
direction of n and negative in the opposite direction. Now at the time to + dt, the location of the surface /(x, to -/-dt) = 0 is 
obtained by moving each point (u 1 , u 2 , 0) on /(x, to) - 0 to {u 1 , u 2 , v n St). Transfer coordinates (u 1 , u 2 ) from /(x, t 0 ) = 0 
to / (x, to + dt) = 0 by this process and retain the definition of coordinate system u 3 as defined for surface / (x, to) = 0 . 
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2- Assume that the limits of the Gaussian coordinates ( u 1 , u 2 ) in defining the moving surface / = 0 is independent of time 
for te[to, t 0 + dt]. We call this two-dimensional domain Q u . See Farassat and Myers ^ for an example. This assumption 
allows us to exchange the order of integration after the second equality sign in eq.(25) blow. 

We can now do all the following manipulations rigorously: 


rto+o t r ( dd) \ rk)+ot r> ^ 2 , 

I I \(f> t yr \ n E\dS dt - - j I v n [3 + vj;) 0 n (u 1 , u 2 , u 3 , t) J 2 ) ( u 1 , u 2 , u 3 , t) du 1 du 2 dt = 

Jto JgnA ° n ' Jdn u 

r rto 

-I I V n (3 + v 2 n ) n ^ g {2) dt du 1 
dd£l„ 


du 2 = 


(25) 


f f Vn (3 + v 2 n )F n ^- dul du 2 =-f f (5 + vJ) 

d d£l u do d Q£l u do 


6 (j) 

dn 


g( 2 ) dn du 1 du 2 = 0 


Here g^) is the determinant of the coefficients of the first fundamental form^ . i n the last step, we have used 
d 0 = 0 (u 1 , u 2 ,v n dt, t 0 +6t)-(j (u 1 , u 2 , 0, t 0 ) - 0 - 0 = 0 by the BC. We can now repeat the process again for all time 
te [0, T] to obtain: 


/ 

dn t (T) 


Edx- 0 


E =-ti + - m 2 =o 


0 (x, t) = constant = 0 (x, 0) = 0 


(26) 


Thus the uniqueness of the solution of the wave equation with Dirichlet BC specified on a deformable body in motion is 
proved. 


3.2.3. The Uniqueness Theorem for an Arbitrary Deformable Surface in Motion With Neumann Boundary 
Condition 


The basic idea of proof remains the same as in the case of the Dirichlet BC with some minor differences. Note that we have 
the Neumann BC 0 W =0 on d£l t . We can show that, after using the BC, we have 


[$ t “ v « - -v «[^ - 2v t (j) T (j) t + (l + v^)0|| 


(27) 


where v- V r 0 = v T 0r where vt and (j)j are the tangential velocity of the surface and the directional derivative of 0 in the 
direction of the tangential velocity, respectively. We now repeat the procedure for Dirichlet BC except that we should start 
from t - 0 and use the fact that 0 W = 0 to show that 0 = 0 is transferred at each step from /(x, t) = 0 to the surface 
/(x, t + dt) = 0, te[ 0, T], so that the right side of eq.(27) always vanishes. Equation (26) and the uniqueness result for the 
Neumann BC then follow. 

Remarks - 1- Since our reasoning is local, the above method can also establish the uniqueness of the solution of the wave 
equation with the Robin BC (a linear combination of the unknown and its normal derivative). 

2- To the best of our knowledge, the above theorems are proved for the first time here. We have extensively searched mathe- 
matics and physics literature and have not found the proofs of these theorems elsewhere. Note that the complications and 
difficulties in the proof come from the motion of the surface. 

3- The above uniqueness theorems ensure the existence of EGF for the problems under consideration. The reason is as 
follows. For a linear operator, the existence of Green’s function is guaranteed if zero is not an eigenvalue of the operator. 
When zero is an eigenvalue of the operator, it means that uniqueness breaks down. Since we have established uniqueness, 
zero cannot be an eigenvalue of the problem and the EGF exists. 
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3.3. The solution of the Lighthill Equation with the Dirichlet or the Neumann BC using 
EGF 


For completeness and for the purpose of comparison of noise prediction formulas in AFN prediction methodologies, we give 
the solution of Lighthill’ s equation using EGF. We assume that the moving surface / = 0 is deformable. The two initial- 
boundary value problems that we will consider are: 

□V = irir- x ^nte[o,n 

OXi OXj 

dp'( x, 0) (28) 

IC: p'(x, 0) = g(x), — — =h(x) xe£l t , 

dt 

BC : p' (x, t) = k(x, t) x e d£l t , t e [0, T] Dirichlet BC 
and 


d 2 T 

D 2 p' = - — xefl t , re[0, T], 

OXi OXj 

dp' (x, 0) 

IC: p'(x, 0) = g(x), — =/z(x) xeFL t , 

ot 

dp' (x, t) 

BC : — — = k (x, t) xedkl t , te[0, T] Neumann BC 

dn 


(29) 


We are considering initial-boundary value problems here. The EGFs for the above wave equations satisfy the following 
conditions: 

U 2 xt) G(x, y 9 t-r) = S(x - y)5(t-r) x ekl t , y e£l t , r < t, t-re[0,T] 

G(x,y,t-r) = 0 xe(l h y e£l t , r > t (30) 

G(x,y,t-T) = 0 xed£l t , ye£l t , t-re[0, T] For Dirichlet BC 

and 

U 2 m G(x, y, t-r) = 6(x - y)5(t-r) x eFl t , y eFl t , r < t, t-re[0,T] 

G(x, y, t -r) = 0 xe(l h y e£l t , r > t 
d G(x, y , t — t) 

= 0 xed£l t , ye£l t , t-re[0,T] For Neumann BC 

dn x 


The solutions of eqs. (28) and (29) are: 


P'(x, t) 

f ( 

Jnxo) L 


fj 

Jo Jn, 


d 2 T u 


g(y ) 


w dy t dyj 
dG(x , y, t-r) 


G(x, y, t - r) dydr - 


P'(x, t ) : 


ff 

Jo Jn t 

f Lj) 

Jn t (0) L 


8t 


d 2 T u 


- h(y) G(x, y,t-T ) 


dy + 


w dytdyj 
dG(x , y, t-r) 
dr 


G(x, y, t - r) dydr - 


- h(y) G(x, y,t-r) 


dy + 


ff 

Jo Jon t 


fl 

Jo Jdftt 


dG(x , y, t-r) 
k(y, t) — dS y dr For Dirichlet BC 

(r) dn y 


(32) 


(33) 


k(y , r) G(x, y, t - r) dS y dr For Neumann BC 


(T) 


In these equations the region fl t (r) is external to the surface dft t (r ) : f(y , r) = 0. We have left out all the details of the 
derivation which involve some subtleties. For example, one has to use the transport identity of eq. (15) in the derivation of 
these solutions. In addition, some of the properties of the Green’s function on either side of r-t must be exploited (See the 
books by Goldstein^ , Stakgold^ and Morse and Feshbach^). 
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We mention here that in acoustic applications, it is usually assumed that the noise generation process startes at t - -oo, or at 
a finite time, after any transient radiation has subsided. Therefore, one rarely sees the volume integrals on the right of 
eqs.(32) and (33), which account for initial conditions (IC), in current noise prediction formulas. This assumption is made 
here as well. 


4. Implications of the Uniqueness Theorem 


From eq.(33), we may conclude that, in the case of a stationary airfoil, when the Neumann BC dp' /dn = 0 is appropriate, 
the acoustic pressure for a TE noise prediction is given by 


p'(x, t ) = 



d 1 T ij 

- — - — G(x, y,t-T)dydT 
dyt dyj 


(34) 


This equation is identical Fighthill’s equation for jet noise prediction, except that Green’s function in the above equation is 
the EGF that satisfies eq.(31). Therefore, as suggested by Howe, we are justified in using incompressible turbulent flow 
simulation near the TE to model the source in the integral above. In addition, statistical analysis on this equation to relate the 
statistics of the noise to the statistics of the turbulence is relatively easy^ . A word of caution is in order here, as the simplic- 
ity of this equation should be balanced against the difficulty of obtaining the EGF, both analytically and numerically. 

Experimental and computational evidence indicate that even at high frequencies the wavefronts of the acoustic waves 
spreading from the TE area are normal to the airfoil surface (see fig. 5) in many TE noise problems. This means that the 
assumption of the vanishing of the normal derivative of the surface pressure is a reasonable one and use eq.(34) is appropri- 
ate for noise prediction at high frequencies and when the airfoil is in motion. Therefore, the use of EGF with the Neumann 
BC dp' /dn-0 offer clear advantages over the use of FSGF in FW-H equation. 

Clearly, eq.(34) can be used to calculate broadband noise from turbulent flows over a thin surface moving in its plane. This 
application has not been explored although there may be situations where such a model can be utilized. 

In the general cases where we cannot use d p' / dn - 0 , we must specify dp' / dn on the surface for the Neumann BC. The 
condition on the EGF is then dG/dn -0 on the boundary of the moving surface. One then does not require the surface 
pressure data as input. But the noise radiation formula, eq.(33) will have both volume and surface integrals just as the solu- 
tion of FW-H equation! There does not appear any advantage of EGF over FSGF. 

When we use the Dirichlet BC, the surface pressure should be specified on the moving surface and the EGF satisfies G = 0 
on this surface. The radiation formula is eq.(32) which has both surface and volume integrals. In this case, the input data for 
both the solution of FW-H equation and eq.(32) are identical. There is the possibility that when EGF is utilized, the cancella- 
tion between the volume and surface integrals is less than in the case of using FSGF^. if such a possibility were proved 
valid, this would imply less sensitivity in an EGF formulation to the accuracy of the input flow data. 

The current experience with the EGF approach has already encountered some limitations in the application of this method in 
aeroacoustics. Doing analytical work on the EGF, and evaluating it numerically on a computer, incurs many problems, such 
as singularities on the surface. Furthermore, for complicated surfaces of interest in AFN prediction, the EGF can only be 
obtained numerically. Therefore, any realistic AFN calculation with an EGF formulation will involve intensive computer 
usage. It should be mentioned here that these calculations are far less than that for turbulent flow simulations (FES and 
DNS). Thus, in general, the use of EGF is within our reach today. Therefore, we think that the use of EGF should be 
explored further for Direichlet and Neumann BCs. 
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5. Exact Green’s Functions for a Wing in Uniform Rectilinear Motion 

Much efforts has been spent to find either analytically^ or numerically^ - ^ the EGF for nontrivial geometries. The most 
significant lesson that one learns from these works is that constructing EGF either analytically or numerically is difficult. 
With the exception of the article by Hu, Guo and Jones^ EGFs have been constructed for two-dimensional geometries. 
However, Dunn and Tinetti have developed the Fast Scattering Code (FSC) to calculate the scattering of the engine noise 
from thin bodies in motion^>51 The FSC solves the Helmholtz equation for a surface in uniform rectilinear motion on 
which the condition dp' / dn - 0 or an absorbing BC is imposed. In a typical FSC simulation, an incident acoustic source is 
specified off the surface to model the propulsion system noise. The scattered acoustic field is calculated by the efficient 
equivalent source method. The problem is solved in the frame moving with the surface so that the complex pressure ampli- 
tude on and off the moving surface can be found in the moving frame when a source with the complex time dependence of 
e l0Jt is specified. Therefore, if the incident noise source is taken as a constant frequency monopole source with unit strength 
fixed in the moving frame, the pressure distribution on the moving surface gives the complex EGF in the moving frame. 

We present results for two three dimensional wings with realistic geometries. The designs and dimensions of these wings are 
shown in figs. 6 and 7. Figures 8 and 9 show the locations of sources for which EGF on the wing surface is computed. The 
radius of the circle of the sources which is located in a plane normal to the wing is the mean chord length for each wing. Each 
source is spaced 10 degrees apart on the circle. The zero degree is the location ahead of the FE, 90 degrees is exactly above 
the suction side and 270 degrees is exactly below the pressure side of the wing. We only present results for three source 
positions for each wing at two frequencies of 3000 Hz and 4097 Hz in the figs. 10 to 13. 

The values plotted in these surface distributions are the absolute value of the complex amplitude of EGF on the wings for the 
observer at the assumed source position. In other words, we are presenting in color 

| G (x, y, t - r) | = | G (x, y) e l w (?_r) | = | G (x, y) | for the observer x at the source position and the field point y varying over 
all the points on the wing surface. Therefore, we have high sensitivity of the noise to the sources on the regions of the wing 
that have high magnitudes (white or red) at the observer position. In particular, we see that for some observer positions, the 
noise is sensitive to the sources over the trailing edge. Note that only a limited range of TE sources contributes to the noise at 
the observer positions shown. This graphical representation of EGF on the wing can serve in identifying the directivity of 
sources on specific areas of the wing by varying x, or given the observer position, we can identify regions of the wing that 
radiate toward the observer at x . 

The FSC is a highly efficient method of calculating EGF of the wave equation with the Neumann BC. It gives aeroacousti- 
cians yet another tool to predict AFN. 


6. Concluding Remarks 

In this paper, preliminary results are presented for a study involving a critical examination of the potential use of EGF 
formulations for the prediction of broadband noise due to turbulent flow over an airfoil surface and its TE. In the process of 
studying this problem we have proved that for a moving deformable surface, infinitely thin or with thickness, the initial- 
boundary value problem of the wave equation has a unique solution for either the Dirichlet of the Neumann BC on the 
surface. This existence of EGF then follows from these theorems. 

Our conclusions concerning the use of EGF are: 

1- For a stationary surface, the use of EGF for TE noise prediction with the Neumann BC d p' / dn - 0 results in a volume 
integral involving the quadrupoles in the vicinity of TE. An incompressible turbulent flow simulation near the TE is suffi- 
cient for low frequencies but compressible turbulent flow simulation may be needed for high frequencies. 

2- It appears that the same Neumann BC dp' / dn - 0 can be used for TE noise prediction for an airfoil in motion at higher 
Mach numbers. The choice of compressibility of turbulent flow simulation near the TE must be decided by the frequency 
range of TE noise and the speed of the airfoil. 
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3- The use of EGF with Neumann BC makes the statistical analysis of TE noise much simpler than when FSGF is used. This 
conclusion applies only when we can use the BC dp' / dn - 0 on the airfoil surface. Otherwise, no advantage over the use of 
FSGF can be realized. 

4- In the situations where the condition dp' I dn = 0 does not apply, the radiation formulas for the Dirichlet and the Neumann 
BC have both volume and surface integrals. No clear advantage of the use of EGF over FSGF can be seen as to the amount of 
input data and ease of statistical analysis. 

5- The use of EGF with the Dirichlet BC has not been fully explored. The input data needed are identical to when FSGF is 
used to solve FW-H equation. The only possible advantage that one can foresee is the reduction of cancellation between the 
volume and surface integrals as compared to the use of FW-H equation. Because of the availability of the method based on 
FW - H pds and the present difficulty of getting EGF, perhaps this advantage will not give enough incentive to researchers to 
examine such an approach. 

6- We have not been able to discern any advantages for EGF for predicting the noise from turbulent BE over an airfoil 
although for a stationary surface, the ease of statistical analysis of the broadband noise when the Neumann BC dp' /dn-0 
is used may be an advantage over the use of FSGF. 

7- We feel that the use of EGF in AFN prediction is still an open problem. We, therefore, recommend that the use of EGF be 
explored further in acoustics. 
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Figure 1. Description of a surface in motion by f{x,t) - 0 , n - V/ , and / (x,t) > 0 outside the surface. 



(OBS.) 




Figure 2. Piston in the wall problem. The darker shaded region is vibrating 




Figure 3. Deformable surface in motion described by f(x,t) = 0 , n = V f , and / (x,t) > 0 outside the surface. 



Figure 4. Gaussian coordinate system on a deformable body in motion. 




Figure 5. Radiation of waves from a TE, showing the wavefronts propagating normal to the airfoil surface. 
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Figure 6. The planform of model 777 wing used for Green’s function 
calculations. The plane of the source (observer) positions is 
also shown. 
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Figure 7. The planform of model Trap Wing used for Green’s function 
calculations. The plane of the source (observer) positions is 
also shown. 


777 wuing scaled lo trap wing span 
local ions of acoustic sources 



Figure 8. The source (Observer) positions around a circle in a plane for 777 
wing. Only the EGF for three observer positions are reported here. 



Figure 9. The source (Observer) positions around a circle in a plane for the Trap 
Wing. Only the EGF for three observer positions are reported here. 
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(a) Monopole acoustic source at 190° from horizontal. 



(c) Monopole acoustic source at 210° from horizontal. 


Figure 10 - Magnitude of acoustic pressure (Pa) on scaled Boeing 777 wing (b/2 = 2.15 m), 

f = 3000 Hz, M = 0.2. 
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(a) Monopole acoustic source at 190° from horizontal. 
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(b) Monopole acoustic source at 200° from horizontal. 
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(c) Monopole acoustic source at 210° from horizontal. 


Figure 11- Magnitude of acoustic pressure (Pa) on scaled Boeing 777 wing (b/2 = 2.15 m), 

f = 4097 Hz, M = 0.2. 
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190° from horizontal. 


Suction Side 




Figure 12- Magnitude of acoustic pressure (Pa) on trap wing high lift system, f = 3000 Hz, M = 0.2. 
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190° from horizontal. 


Suction Side 




(c) Monopole acoustic 



source at 210° from horizontal. 
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Figure 13- Magnitude of acoustic pressure (Pa) on trap wing high lift system, f = 4097 Hz, M = 0.2. 


